General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



v_/ I 


EY6- 1 03 05 

/ i G? ' lZ^ 

^ » 5 i . vO T W i. C L\ \,'V^ 'Wv ("T ' <r 

' ' I V. . 

*’?-'j'f3‘ J^siht'e ci'.Il, i</-.. . ■ . ■ . 

in too interest of early aru; /::•■ •.':. 
semination tf Earth Resourcfs 
Program Information and f,. ji! H 
fir any uSe made Ihereoi.’' 



I^.SP'C-CVo f~{ f 


5f 


CN 


VO 


CN 

U) 

CN 

(tj in 

• 

■H o 

•o 

U ro 


c o 


3 O 


ro 

\ 

ro 

U 

a sj 
nJ CN 

0) o 
• a< 

U O K-l 
0) M O 
• • w p VI 
CN W U 



E-! 4J Q) 
tx, u ja 
•a! O -P 
Oj 0< 

0) *W 
« O 


QUATERLY REPORT 



E >t C 
«« iH O O 
» M .H O 
CJ 5 0 ) W • 

O -P w m 
oa Li -ri m 
(0 a 
S B (j 
M O U E 
E-i U 
in — Qu 
w to 

m w VO cj 
ca H r~ 'j> 

ptfl M CTi 

to T— 

I » 

^t-t • w 
m to a 0) 

O W rtJ -i-t 
m >-3 p 

0 H 

^ E I C 

1 U 3 
VO E m s 
r' w B 
w P 3 cr o 

W UK T— cj 


DECEMBER 1975 - JANUARY 1976 

^8 79<0 

3 \S| 902.6 



T 





O/lginal photography may be purchased from: 


EROS Data Center 
10th and Dakota Avenue 
Sioux FallSt SD 571SS 


oraemiu. cohtmhs 

BOIOR ILLBSTBATlOtiS 


The Quaterly report of the french part of the AGRESTE 
program presents the directions in which the different laboratories 
will carry out their investigations. It shows the ground-truth works 
at the time of the Landsat 2 passes on the rice-fields of Caraargue and 
the ground-truth methodologic works on poplar-site. The different kinds 
of imagery were available on November 1975 and the first data processing 
in this report corresponds to the preliminary investigations. 

The different chapters show the interest that the rice- 
field temporal mapping presents, the possibility of alluvial facies 
localization using infra-red color photographs, the methodology of 
measuring cover rate in poplar stands with a first inventory and finally 
the specific results concerning every kind of remote sensing on the 
poplar-site. 


G, FLOUZAT 


Technical Manager 
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CHAPTER I 


FIRST RESULTS FROM THE AGRESTE PROGRAMME 


ON IRRIGATED CROPS AT FRENCH TEST-SITE n“ 4 


R. MARIE 

Station d' (melioration des ylantes - INRA-Montpellier 


T. LE TOAN 

Centre d’ Etude Spatiale des Rayonnements 
Toulouse 


INTRODUCTION. 

The objectives of the Agreste program on irrigated cultivation 
(rice) in the Camargue region French test-site n° 4 (figure 1) are the following : 

- to make an inventory of the total area of rice field within the 

test site, 

- to estimate production, 

- to identify plant diseases. 

To date, only the first objective of the program has been studied using LANDSAT 1 
imagery, since the data were insufficient to cover the growth stages of the rice 
during one vegetation cycle. 

Research has been carried out along the following two lines ; 

1. Determining the necessary ground truth interpreting the 
LANDSAT imagery. This work is conducted by the Station d ’Amelioration des 
Plantes (Plant- improvement station) of the INRA in Montpellier. 

2. Traitment of the LANDSAT imagery, carried out by the Centre 

d 'Etude Spatiale des Rayonnements in Toulouse (Space radiation study center). 

In this paper these two aspects will be discussed and the first 
results presented. 




FIGURE 1 


FRENCH AGRESTE TEST SITE N° 4 



1 - GROUND DATA DETERMINATION. 


Information was collected on the state of rice fields and on climatic 
factors in the surrounding areas during the passages of the LANDSAT 2 satellite 
for the vegetation season of 1975. 

The state of the rice fields is determined by ground observations 
which involve the stages of development of the rice for the three varieties 
cultivated in one parcel knoxm as the Mas de Gouine (figure 2). These are 
the three most commonly cultivated rice varieties in the Camargue region, 
i.e.. Delta, Euribe, and Balilla 28. 

During the passages of LANDSAT 2, observations on the stages of 
development and on soil conditions were made (Table 1). The legends to the 
abbreviations corresponding to the stages of development are as follows ; 

LE ; (levee) emergence 

DT ; (debut de Tallage) : Beginning of the Tillering 

DM ; (debut de montee) : Beginning of the Stem Elongation 

GO ; (gonflement visible) : Boots just visible 

DE ; (debut d’epiaison) : Beginning of Heading 

ME : (mi floraison) : Mid flowering 

LA : (grain laiteux) : Early milk stage 

PA : (grain pateux) : Early dough stage 

DU : (grain dur) : Caryopsis hard 

FC ; (fin de cycle) Caryopsis loosen, ing in daytime. 

At the same time, climatic measurements were made during the pas- 
sages of LANDSAT 2 by meteorological recorders placed in an experimental rice 
field (the Mas d ’Adrien, Figure 3) with the aid of the. Service Technique 
d ’Etude des Facteurs Climatiques de 1 ’Environnement (SIEFCE-INRA, Montf.ivet) 
(the Environmental Climatic Factors Technical Study Service). 

The measurements given in Table 2 are as follows : 

1. Air temperature : recorded at 2 meters in meteorological shelter 
conforming to international standards. Three values are recorded over a 24- 
hour period : maximum, minimum, and average. 
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FIGURE 2 

"MAS BE GOUINE"LOCATING 
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FIGURE 3 

"MAS D ’ADRIEN" EXPERBIENTAL PARCEL 





n“ 

Date of passage 

DEVELOPMENTAL STAGES 

Soil conditions 

DELTA 

EURIBE 

BALILLA 28 

n 

2 MARCH 

•• 



DRY 

2 

20 MARCH 

- 

- 

- 

WET 

3 

7 APRIL 

- 

- 

- 

WET 

4 

25 APRIL 

- 

- 

- 

DRY 

5 

13 MAI 

LE 

LE 

LE 

UNDER WATER 

6 

31 MAI 

LE 

LE 

LE 

UNDER WATER 

7 

18 JUNE 

DT 

LE 

LE 

UNDER WATER 

8 

6 JULY 

DM 

DT 

DT 

UNDER WATER 

9 

24 JULY 

GO 


DM 

UNDER WATER 

10 

11 AUG. 

MF 

DE 

GO 

UNDER WATER 

11 

29 AUG. 

LA 


MF 

UNDER WATER 

12 

16 SEPT. 

PA 


LA 

UNDER WATER 

13 

4 OCT. 

FC 


PA 

PARTLY AT ZERO LEVEL 

14 

22 OCT. 


FC 

DU 

- 

15 

9 NOV. 



ED 





TABLE 1 : DEVELOPMENTAL STAGES OF RICE PLANT OBSERVED AT EXPERIMENTAL 


FIELD 














AIR TEMPERATURES 


DATE OF 
PASSAGE I Mini Maxi 


75.03.02 

75.03.20 

75.04.07 

75.04.25 

75.05.13 

75.05.31 

75.06.18 

75.07.06 

75.07.24 

75.08.11 

75.08.29 

75.09,16 

75.10.04 

75.10.22 

75.11.09 



13.5 
8,5 

15.5 

29.0 

21.7 

17.0 
21,9 

25.0 

30.7 

28.5 
25,2 

24.8 
27,7 

18.0 
9,4 


12,5 

4,5 

10,3 

19,1 



15.1 

15.9 

20.9 
24,0 
23,3 
20,8 

19.2 

20,6 

12.3 
7,5 


8,9 

4,2 

14,2 

15.8 
15,6 

14.0 

12.0 

11.9 
3,0 
2,8 


15.8 

16.9 
16,4 
16,0 
13,0 

12.9 
5,0 
3,7 


SOIL TEMPERATURES 

a-lC 

cm 

a-50 

Mini 

Maxi 

cm 

9,5 

10,0 

9,6 

7,3 

8,9 

8,8 

8,0 

9,5 

9,0 

14,5 

15,1 

15,1 

11,4 

15,5 

14,1 

15,8 

19,4 

18,4 

15,4 

20,0 

20,4 

16,9 

20,2 

20,9 

20,7 

26,4 

24,2 

21,4 

26,5 

26,0 

16,4 

22,7 

21,4 

14,0 

19,5 

20,5 

13,0 

20,7 

19,7 

8,0 

13,6 

13,5 

6,0 

10,2 

12,5 


SUN RELATIVE 

(h) HUMIDITY 
<40 >80 
% % 



0,0 

49,8 

0,2 

0,02 

0,0 

11,1 

0,0 

0,0 

0,0 

0,0 

0,1 

0,0 

0,0 

0,2 

18,3 


TABLE 2 : CLIMATIC MEASUREMENTS RECORDED AT THE "MAS DE GOUINE" 


AREA 


























2. Minimum actinothermal index : recorded by thermometers at 

10 and 50 cm below the soil. 

3. Wind ; prevailing direction ( 1 ; HE, 2 ; E, 3 : SE, etc 
and average velocity at 2 meters above the ground, measured by anemometer 
(expressed in meters per second) . 

4. Soil temperature : measured at - 10 and - 50 cm beneath short 

grass . 

5. Sun : expressed in hours and tens of an hour using Campbell 
heliograph recordings. 

6. Relative humidity , Expressed in j 

a. Number of hours during which the relative humidity was less 
than 40 %, 

b. Number of hours during which the relative humidity was over BO %. 

7. Rainfall : expressed in mm per day, obtained by a rain gauge. 

In the meantime, experiment is carried out for’-studying the growth 
stages of several rice-plant varieties : Balilla 28, Cigalon, Cristal, Delta, 
Euribe when the sowing dates are intentionally spread out in time. 

The earliest date was April 21 and the latest, May 05. Table 3 gathers 
the evolution of each variety allowing to one of the 4 sowing dates. On can 
note that the difference of growth stageSjimportant in the first phases, 
decreases more and more to harvest. Further investigations will be done to 
study the relation between climatic conditions and the growth stages with the 
data from Table 2 and Table 3. 

11 - PLANNING OF LANDSAT 2 DATA PROCESSING, 

In conjunction with ground observations and measurements, a plan of 
the treatment of the LANDSAT data can be made for each passage of the satellite. 

The first stage of such treatment involves rice field mapping^rice- 
varieties mapping and diseases detection (Table 4) . 


VARIETY 

SOWING 

First 

leaf 

BALILLA 28 


07.05 

CIGALON 


07.05 

CRYSTAL 

21.04 

07.05 

DELTA 


07.05 

EURIBE 


07.05 

BALILLA 28 


12.05 

CIGALON 


12.05 

CRISTAL 

29.04 

12.05 

DELTA 


12.05 

EURIBE 


12.05 

BALILLA 28 


15.05 

CIGALON 


15.05 

CRISTAL 

07.05 

15.05 

DELTA 


15.05 

EURIBE 


15.05 



08.06 

08.06 



04.06 


06.06 
06.06 
06.06 
06.06 
06 .'06 


Beg. of 1 
Stem 
Elon- 
gation 

Boots 

just 

visible 

25.07 

06.08 

12.07 

23.07 

25.07 

06.08 

20.07 

26.07 

22.07 

06.08 

20.07 

31.07 

12.07 

20.07 

20.07 

29.07 

12.07 

20.07 

26.07 

02.08 

28.07 

03.08 

12.07 

19.07 

26.07 

02.08 

15.07 

22.07 

30.07 

06.08 


BALILLA 28 


25,05 

10.06 

28.07 

04.08 

CIGALON 


25.05 

10.06 

12.07 

19.07 

CRISTAL 

15.05 

25.05 

10.06 

30.07 

04.08 

DELTA 


25.05 

10.06 

17.07 

25.07 

EURIBE 


25.05 

10.06 

25.07 

03.08 


TABLE 3 : DEVELOPMENTAL STAGES OF 5 RICE PLANT 


Heading 

to 

flowering 


End 

of 


Early 

milk 


flowering stage 


12.08 

31.07 

10.08 

03.08 

11.08 


09.08 

30.07 

06.08 
02.08 
11.08 


09.08 

31.07 

09.08 

03.08 

12.08 


09.08 

31.07 

08.08 

05.08 

12.08 


20.08 

08.08 

20.08 

09.08 

23.08 


14.08 

08.08 
12.08 

09.08 

19.08 


18.08 

08.08 

16.08 

11.08 

21.08 


21.08 

09.08 

16.08 
12.08 
23.08 


27.08 

13.08 

25.08 

20.08 
04.09 


27.08 

13.08 

18.08 
L8.08 
29.08 


25.08 

13.08 

23.08 

21.08 
28.08 


27.08 

18.08 

24.08 

21.08 
04.09 


Early 

dough 

stage 

04.09 
21.08 

06.09 

02.09 

15.09 


03.09 

21.08 

27.08 

25.08 

09.09 


09.09 

21.08 

30.08 

27.08 

09.09 


11.09 

25.08 

03.09 
28.08 

18.09 


Good 

|Caryopsis for 
hard I harves t 


18.09 

08.09 


08.10 

22.09 


21.09 

15.09 

28.09 


03.10 

27.09 

05.10 


OBSERVATIONS 


Ahnormaly long 
cycle 

Heterogeneous 

emergence 


22.09 

07.09 

15.09 

09.09 

25.09 


05.10 

22.09 

30.09 

22.09 

05.10 


25.09 

09.09 

20.09 

11.09 

22.09 


14.10 

22.09 

03.10 

25.09 

07.10 


25.09 

12.09 

16.09 

10.09 

30.09 


14.10 

25.09 

06.10 

25.09 

15.10 


Disease 


VARIERIES AT STAGGERING SOWING DATES 



















DATES 

OBJECTIVES 

MAY 13 

Rice-field Mapping 

MAY 31 

Rice-field Mapping 

JUNE 18 

Rice-field Mapping 

JULY 6 

Rice-field Mapping 

JULY 24 

Rice-field Mapping 

AUGUST 11 

Rice-field Mapping 
Disease detection 

SEPTEMBER 16 

Rice-field Mapping 
Disease detection 

OCTOBER 4 

Mapping of rice varieties 

OCTOBER 22 

Mapping of rice varieties 


COMMENTS 


Results including all natural bodies 
of water. 

Difficulty : rain gauge reading 11 mni 
Difficulty : a certain number of rice 
fields were dried out for weed-killer 
treatment. 

Results including rice-fields and 
marshes 

Results including rice fields, marshed 
and temporarily flooded plots (vineyards 
wheat-field after harvest ...) 


TABLE 4 : PLANNING OF LANDSAT 2 DATA PROCESSING 
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1. Mapping of surfaces under water : determined by the classifi~ 
cation of imagery on the following dates s May 13, May 31, June 18, during 
which the rice fields VJe^re under water. This poses several problems : 

- Natural bodies of water give the same response as rice fields, 

- May 31 is distorted by the appearance of rain (rain gauge reading : 

11 mm) , 

“ One June 18, certain areas were dried out for various chemical 
weed-killer treatments. 

2. Mapping out rice fields in the vegetation stage : the July 6, 

July 24, and August 11 passages can be used to identify the rice under 
following states ; emergence to tillering, stem elongation , booting ^appea- 
rance of first leaves above the water ; in addition, rice fields can be 
distinguished from other areas of vegetation in the MMS7 band, and from natural 
bodies of water in the MSS5 band. 

Certain difficulties may arise : 

The marshes, which are extensive, particularly in the southern part 
of the zone under study, can appear as rice fields (grass on marsh). 

On dates at the end of July (July 24) certain fields were purposely 
flooded by users (vinegards, orchards). 

3 . Mapping out rice varieties : 

The October 4 and 22 imagery can be used to discriminate between 
different varieties of rice. 

On October 4, the Delta variety was at the end of the vegetation 
cycle and the plants were yellow ;the Euribe variety was starting to turn 
yellow; the Ballila 28 was still green. 

On October 22, most of the rice fields of the Delta variety had 
been harvested; the plants of the Euribe variety were yellow; the Bailila 28 
pla..ts were just beginning to turn yellow. 

4. Detection of diseases ; 

The passages on August 1^ and 29, October 3 and 4, can be used to 
detecf diseases which can be found at different points in 1975 , 
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The disease is caused by PYE.ICULAE.IA ORYZA on parcels over treatment 
with nitrogen fertilize^ when occur particularly wet and harm climatic conditions 
within several days. 

The visible signs of the disease are detected on the field from 

August 10 

III - LIST OF LANDSAT 2 IMAGERY AVAILABLE OVER FRENCH TEST-SITE n” 4. 


The LANDSAT 2 images received today are the following: 


n° 

DATE 

IDENTIFICATION 

OBSERVATION 

1 

JUNE 17 

2146-09453 

Cloud covering upon test-site - Unusable 

2 

JULY 4 

2163-09394 

Unusable 

3 

JULY 5 

2164-09453 

Unusable 

4 

JULY 6 

2165-09511 

Good - Compatible - Tape ordered 

5 

JULY 23 

2182-09451 

Good 

6 

JULY 24 

2183-09505 

Good - Compatible Tape ordered 

7 

AUGUST 11 

2201-09502 

Good - Compatible Tape ordered 

8 

AUGUST 29 

2219-09500 

Unusable 

9 

OCTOBER 3 

2254-09440 

Good - Compatible Tape ordered 

10 

OCTOBER 4 

2255-09495 

Good » 


On the whole, 6 imageries are available 2 of which are redondant 
(July 23 or 24, October 3 or 4) . 

Allotting for these available data, processing is planned for ; 

1. Mapping of rice fields (july 6 and july 23 scenes) 
by multispectral data classification 

2. Mapping of rice-fields and disease detection (august 11 scene) 
by multispectral data classification 

3. Mapping of rice-field by raultitemporal data classification (july 6 
and 23, august 11) to eliminate the errors caused by the complexity 
of particularities: presence of marshes, existence of intentional 
flooded fields other than rice fields, and the case of rice-fields 
dried out for chemical treatments. 

4. Mapping of rice field into subclasses of different rice varieties by 
classification of data on October 3 scene. 
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IV - DATA PROCESSING. 

The data processing developed by the Centre d*Etudes Spatiales 
des Rayonnements can be devided into 5 main parts : 

- data management 

- image display 

- sampling 

“ classification 

according to the chart in figure 4. 

t 

1. Data management . 

The data come from either magnetic compatible tapes, or from digi- 
talization of black and white transparencies on the Joyce-Loebl SCANDIG 25 
microdensitometer. 

In this case, each image is canned row by row with a step of 25 p 
giving a numerical value from 0 to 255 par byte for each elementary surface 
25 p X 25 p. Each LANDSAT imagery (MSS4 to 7) corresponds to 4 files of 
2500 X 2500 = 6,2 X 10^ bytes. 

Many problems result from the vast amount of data, and the solution 
may satisfy these conditions: 

- reduction of the access to information time 

- optimal data handling in a machine word (CDC 7600) 

- perfect coincidence between the digitalized geographical zone 
in 4 MSS transparencies. 

In this order image is devided into elementary zones of 60 x 60 points. 

2. Display. 

For the elementary restitution of the image, optical density histogram 
is studied by these methods : 

- density slicing 

- enhancement by equal probability quantisation and estimation of 
entropy 

- extraction of normal elementary distribution. 
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FIGURE 4 DATA PROCESSING 










3. Sampling : 


For supervised classification methods selected samples must be located 
according to ground truth determination and aerial photograph interpretation. 

The process is refined to improve the homogeneity of the sample 
defined by their statistical characteristics. 

4. Data reduction : 

Necessary in the case of parameters redundancy, or when only 
certain parameters are specifically adapted to the type of classification. 

5. Classification : 

The classification model identifies a ground pattern resolution 
element represented by the feature measurement vector. Several models are 
used : 


1. Minimum distance to the mean optical density measurement 

2. Elliptical boundary criterion 

3. Maximum likelihood ratio 

4. Clustering method. 

V - PRELIMINARY RESULTS FROM LANDSAT I IMAGERY. 


The LANDSAT I images used for the preliminary investigations in 
the AGRESTE program concern 2 dates : March 21 and May 14, 1973. The purpose 
of the treatments applied is to test the proposed methods. The data do not 
correspond to the vegetation period of the rice and the aim of the data 
processing is to detected damp ploughed parcels in the month of March and 
flooded parcels in the month of May. 

The first mapping was carried out using the May 14, 1973 imagery, 
when the rice fields were under water. (The flooding was announced in all 
of the parcels except for 1 or 2 in which the flooding took place May 15) . 


Samples are taken on areas located by observations on the ground 
determining rice parcels, vineyards, cereals, as well as some representative 
samples of the class "Rhone" and "cloud". The study of the responses to these 
samples in the MSS 5 and MSS 7 bands (figure 5) reveals the grouping of the 
following major classes : 

water bodies 

vineyards 

cereals 

In the class of water bodies two sub-classes can be distinguished ; 
rice fields and the Rhone. 

By studying these responses, it seems possible to reveal rice fields 
by means of density slicing method. The result (figure 6) shows that it is not 
possible to distinguish the responses of ricefields from those of the Rhone. 

Thus the result includes all vjater bodies in the zone. It can further be noted 
that the "petit Rhone" in this case has the same density as the vineyard samples. 

The multispectral classification, from MSS 4 to 7 , using supervised 
methods, is then applied to the zone. Figures 7 and 8 show the results of the 
mapping out of a zone 180 by 180 elementary points (an area of about 12.5 
by 12.5 km) . 

Several things can be noted from the output listing : 

“ an East-West distortion due to the characters of the printer, 

- each elementary spacing corresponds to a ground surface of 80 x 80 m 

- similarity between graphic symbols chosen for the Rhone and rice 
fields, 

- clouds are left blank. 

On the result in figure 7, showing a part of the test-zone, the 
following can be noted : 

- both the stream (petit Rhone), and the edges of the river (Rhone) 
are classed with the rice fields. The first conclusion to be drawn is that on 
that date, May 14, rice fields cannot be distinguished from natural water bodies. 


MSS 7 






FIGURE 6 

DISPLAY BY DENSITY SLICING 
THE 14 MAY 1973 SCENE 





FIGURE 7 



CLASSIFICATION OF THE 14 MAY 1973 SCENE 







~ urban areas (the to\ra of Arles) have similar responses to those of 

vineyards, 

“ the test consisting of subdividing the class of rice fields is 
inconclusive. 

In figure 8 the vineyard and urban classes have been removed from 
the results. The same symbol is given for rice fields, the Rhone and the 
Petit Rhone. 

The second mapping is obtained using the MSS 7 imagery of March 21, 
1973 and May 14, 1974, The purpose of this multitemporal classification is as 
follows : 


- to test the coincidence between images taken on different dates, 

~ to eliminate confusion between rice field and water body responses, 

- to class points below the clouds in one image or the other. 

Figure 9 shows the study of samples taken at two different dates. 

It can be noted that the classes studied, i.e., rice fields, the 
Rhone, the petit Rhone, urban areas, vineyards, cereals are separated on the 
figure. 


The result of the classification in figure 1 0 shows that : 

- it is difficult to distinguish between rice field responses and 
those of the petit Rhone and the outline of the Rhone, which proves that the 
samples studied in figure 5 do not cover all limiting cases, particularly for 
the outline of the river. 

“ the urban area is still mixed with the vineyard class. 

The conclusion drawn from this result is that the periods of March 
and May are not adequate for discriminating between the above mentioned classes, 
although the goals fixed for this multitemporal mapping have been achieved. 

Figure 11 shows the same result in which the three major classes 
appear : water bodies, vineyards, and cereals or grass. This output was made 
with the Calcomp tracer where water bodies are showed in blue, vineyards in 
green and other surfaces (cereals , grass ...) in red. 


MSS 7 -21 MARS 1973 



FIGURE 9 


KSS7 SAMPLE RESPONSES, 21 MARS 1973 VERSUS 14 MAI 1973 
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FIGURE 10 


MULTITEMPORAL DATA CLASSIFICATION 
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FIGURE 11 

MULTITEMPORAL DATA CLASSIFICATION (CALCOMP TRACER DISPLAY) 







VI - PRELIMINARY RESULTS OF THE AGRESTE PROGRAMME USING LANDSAT 2. 

Studies have been made on the earliest imagery received, i.e., 
from July 23, 1975. 

Figures 12, 13, 14, and 15 show the MSS 4 to 7 images of July 23, 

1975 (2182 - 09451 scene). 

In order to determine samples relating to rice fields on these 
images a study of the aerial documents is carried out, compiled by ground 
truth investigations. 

1. Study of the aerial imagery ; 

The aerial images come from a project organized by the GDTA 
(Groupement pour le developpement de la teledetection aerospatiale) of 
the ONES (Centre National d’ Etudes Spatiales) on June 20, 1975. 

The plane used for the project is the B17 of the IGN (Institut 
Geographique National) flying over the area at an altitude of 8,000 meters. 

The data used within this study consists of images reconstructed 
using 10 “ channel Daedalus Multispectral scanner recordings and specifically 
the two channels closet to the MSS 5 and MSS 7 bands of Landsat : 

- channel 7 ; spectral band from 650-700 nm ; image given in figure 16 

- channel 10 : 9000-1100 nm, figure 17. 

The imagery obtained has a resolution of about 20 m in the central 
part along the axis of flight, for a field of view of the MSS Deadalus of 
2.5 mrd and for a flight altitude of 8,000 m. 

In order to determine the stages of development of the rice at 
June 20, 1975 we refer to the observations made on June 18, 1975 (Table 2), 

On this date, all the parcels are under water and different varieties 
of rice exhibit different stages of development ; Delta is at the beginning 
of the tillering? the other varieties, Euribe and Balilla 28 are at the stage 
of the emergence (first leaf visible). 
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FIGURE 12 


MSS4 LANDSAT 2 FROM 23 JULY 1975 SCENE - SCALE = 1/200.000 











FIGURE 13 

MSS5 LANUSAT 2 FROM 23 JULY 1975 SCENE - SCALE = 1/200.000 
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FIGURE 14 

MSS6 LANDSAT 2 FROM 23 JULY 1975 SCENE - SCALE = 1/200.000 
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FIGURE 15 

MSS7 LANDSAT 2 FROM 23 JULY 1975 SCENE - SCALE = 1/200.000 
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FIGURE 16 


DEADALUS IMAGE FROM 650-700 run BAND - 20 JUNE 1975 
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FIGURE 17 


DEADALUS IMAGE FROM 1000-1100 nm BAND 


20 JUNE 1975 
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In figure 17 , which shows the image taken in channel 10 (near IR) , 
xre note the following ; 

- the optical density of the parcels depends essentially on their 
humidity. Thus the surfaces covered with water are dark on the document. The 
rice fields with Euribe and Balilla 28 are less dark than the Delta rxce 
field, where the presence of more leaves reduce the reflectance of the surface. 

“ within a great area occupied by rice fields, there is a netx^ork 
of canals and dikes. On the satellite images, these areas which do not contain 
rice will be registered as rice fields. This fact will be taken into account 
for estimating the total area of rice planted, using cartography data. 

In addition, on the date of the passage of the plane, several parcels 
of rice fields were temporarily dried out for purposed treatment. 

In figure 16, which corresponds to the image taken in the 650-700 nm 
band (= MSS 5 of LAND SAT) , we note the following ; 

- the areas covered with vegetation are dark on the documeni; ; bare 
land is light. Thus, corn, cereal, soja, and grass fields appear dark, and 
trees along the river are black, 

- the rice fields appear in intermediate grey tones, and the diffe- 
rences in density are observed in the parcels of rice located in figure 17 bis. 
Investigations can be made on specific areas in order to establish correlations 
between such differences in density and differences in variety or differences 
in sowing dates . 

Studies are in progress for automatically classifying the Deadalus 
data in order to perform a precise mapping in a small area of the test zone 
for two purposes : 

- to acquire ground truth for zones outside of the prospection area, 

- to compare the results obtained by both types of data collection, 
i.e., satellite, and airplane. 

2. Sampling: 

The investigations to determine ground truth are in progress. At 
first, such a study is limited to the parcels of rice cultivation only. 

By means of the aerial imagery and ground-based investigation, 
several rice fields are located on a map made by the IGN (figure 18) , These 
samples are then studied on the Landsat imagery for supervised classification. 















3. Treatment of the Landsat 2 data. 


The samples are located on the displays obtained by the density 
slicing of the histogram. Figures 19 and 20 show the displays of one part 
of the area using the MSS 5 and MSS 7 images. 

We notice that in figure 20 the black spots marked 1 on the display 
correspond to an area covered with water. By ground-based investigation, it 
appears that this was a harvested wheatfield which had been flooded after 
burning in order to facilitate ploughing. 

A field identified as a rice field (les Pebrieres , marked 2 in 
figure 20) has a different response in different bands from the other rice 
fields taken for verification. 

These anomalies prove that in order to be identify all rice culti- 
vation areas in the Caraargue, given the complexity of agricultural practices 
in this area, several dates are required for multitemporal classification, 
as well as a thorough knowledge of the ground truth. 

The samples of rice cultivation taken are used to perform an analy- 
tical classification, the purpose of which is to classify rice cultivation 
only. 


On the first rice classification results, showed in figure 21, vje 
note that rice fields identified by ground truth investigations on figure 18 
are well classified. But the results include also other non rice areas which 
behave like rice field at that date (stuff on water) . 

Data processing is in progress concerning multitemporal classifi- 
cation as well as disease detection. 
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FIGURE 19 


MSS5 LANDSAT 2 DISPLAY PROM 23 JULY 1975 SCENE 
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FACIES DISTRIBUTION WITHIN ALLUVIAL DEBOSITS 
INFRARED REMOTE SENSING RESULTS 
B. DONVILLE 


Abstract : 


In fluviatile deposits, the distribution of granulometrecally 
different sediments (clays, sands, gravels) resulting from the flooding 
of former meanders and their development, is very difficult to study. The 
usual methods used are geophysical prospecting or core drills. 

The use of false-color aerial photography is shorn to achieve 
similar results much faster using a new method for observing these docu- 
ments (‘by using a system of Colored lights making it possible to alter 
each spectral response) . 

This paper describes study using this method in the lower plain 
of the Garonne, near St Jory (Haute Garonne) , where our results were 
checked againt local hydrogeological and geophysical data. 

It has thus been possible to reveal a zone of better permeability 
grouping several sandy cordons hurried beneath superficial liraons. 


Laboratoire de Geologic et Geochronologie 
CEARN 

Universite Paul Sabatier 


INTRODUCTION 


The localization of different facies within alluvial deposits has 
a highly practical importance both for hydrogeological problems (water flow, 
drainage) and for primary material problems (gravels ...) This distribution 
is the result of the location of former meanders and their development. Thus 
the first part of this paper will be a brief review of the flooding patterns 

of present streams. 


Generally, corresponding prospection studies are the result of 
specific information obtained either from core drills or from shallow geophy- 
sical prospecting (1). Obviously, in order to obtain sufficiently coherent^ 
data on these highly heterogeneous alluvial formations over a large area fine 
mesh campaigns must be available, Such work is therefore long and costly. 

Near infrared photographs (false color or IRC) can be used very profitably 
for this purpose when the spectral responses are analysed by an adequate 
apparatus. The second part of this article will describle the corresponding 
experimental apparatus after a review of the colorimetric composition of the 

IRC documents. 

In order to evaluate the method, we compare our results with tradi- 
tional processes, taking the lower plain of the Garonne near St Jory (Haute- 
Garonne) as an example. This will appear in the third part of the article. 

Reviev? of alluvial sedimentation 


The formation of alluvial plains depends on the development of the 
balance profile of the body of water and of the extent of its flooding. For 
a certain state of these two components, the outline of the bed will become 
meander-shaped (2). The sedimentary deposits will have varying forms (see^ 
figure lA) depending on the speed of the flow of water in the curve of this 
curvilinear charnel. The concave bank will be strongly affected and will ^ 
thus develop a steep bank. The converse occurs as the distance from the axis 
of the channel increases, fines and finer sediments are deposited on the 

convex bank . 


The meander may extend along its concave side depending on the 
development of its equilibrium profile (figure 1-B) . Thus successive deposits 
would be found in the same order of granulometric classification. Generally 
this movement is uneven so that between each deposit of the convex bank 
there appear shallow depressions giving a crescent shapped structure to 
the whole. When the converse occurs such as when the river tends to resume 
a more rectilinear pattern, the former chenal will be filled with increasingly 
finer sediments, and even more so as the former body will be abandoned more 
rapidly. At the final stage the abandoned body (dead arm) will form a depres- 
sion. 


All such depressions will provide preferred locations for silty 
sedimentation during seasonal flooding. This process as a whole thus provides 
a succession of coarse deposits (sands) , in crescent form, alternating with 
finer (clay) areas (figure 1-C) . 

If this succession is observed longitudinally, one can then notice 
that there will be an average distribution of coarse formations which is less 
windy than the successive patterns of the stream channel, (figure 1-D) . These 
deposits thus formed are practically always, and this is so for the alluvial 
plains of the streams of our region, entirely covered by flood basin deposits 
(clays). The former morphology is thus supposed to be covered up. However, 
due to the differential compaction between clays and sands, the slight alti 
metric variations reveal this former pattern on the surface (figure 1-E) . This 
pattern can be complicated by the superimposing over a period of time, of si 
milar cycles involving a superimposing in any order of the structures previously 
described. 


All such granulometric variations within the alluvial deposits ob- 
viously cause changes in permeability. Hydrogeological research an such deposits 
will therefore consist of delineating the areas having a high concentration of 
sandy meanders (point bar) . These do not necessarily correspond to areas of 
hydraulic drainage related to the form of the substratum. 


-A- 


Treatment of Infrared Color Photographs (I.R.C.) ; 

Aerial photographs in the near infrared region (0.7-1 p) of the 
surface of the Earth, can be obtained either in black and white or in color. 
The method for analysing x-yhich we describe later will show our reasons for 
preferring the latter form. 

These infrared color photographs are called "false color" because 
the document is obtained by means of a spectral shift in the signals emitted 
by the ground. In particular, infrared information xjill be shifted to red and 
the red information to green (3) . 

In figure 2 we have given the effects of this spectral shift for txTO 
sets of information, (red and I.R.) and (green and I.R.). It can be noted, for 
example, than on the photograph a x^hole series of red varying with the vegetal 
density (real green) and x%rith humidity (x'rater absorbs I.R.) will be obtained. 
Thus for a given development of vegetation, we have a means of relating photo- 
graphip information to soil humidity. In addition, it is well knoxra (4) that 
the reflectance of the soil, in this wavelength band, increases as the size 
of its elements decreases and vice versa. Thus sandy soil will appear darker 
than clayey soil for the same degree of humidity. 

If x^e noij illuminate the photograph with a red light and transcribe ■ 
the result into black and x/nite, we note that we xiTill have a darker grey if, 
for example, the I.R. information is low. One can argue analogously for green 
lighting (figure 2, on the right). Thus we have just seen hov7 to transcribe 
certain information obtained very locally on IRC photographs into black and 
white. Now, the layout of cultivated parcels on alluvial deposits provides 
multicolored mosaics on a photograph, making small scale observations very 
difficult. In order to equalize the average of this information on the sur- 
face, it is necessary to be able to amplify or subdue either the red or the 
green information so that dark zones, for example, can be traced from one 
parcel to another. This is obtained in the montage in figure 3. The IRC 
photograph is illuminated by means of two flood lamps, one red and one green, 
each pox^ered by a variable transformer so that the brightness can be adjusted 
individually. The document thus illuminated is televised by means of a camera 
on a tripod and displayed on a monitor screen. This makes it possible to 
analj'se black and X'jhite images on xvhich one can vary the definition (via the 
camera's focus), the brightness and the contrast (the TV screen). 
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Using this device we have treated the 74-2-IRC-257 photograph fur- 
nished by IGN (altitude : 5600 meters, time : 13:10, focal length 152/7) on 
the lower plain of the Garonne west of St Jory. Slightly enlarged, the docu- 
ment, illuminated in red only reveals a dark area lying parallel to the Garonne. 
The boundaries of this area are easier to distinguish when slightly out of 
focus; 


With further enlargement, and with only low green lighting in order 
to avoid over darkening the photograph winding long line patterns are revealed. 
Actually, such observations are better revealed by a varying green brightness 
than by using a fixed value of the variable transformer. 

Finally, by means of our method of green and red ’’’balanced" lighting 
we can observe a series of concentric, primiarily curvilinear darkening zones , 
These results are sketched in figure 4. 

Comparison with ground truth : 

In order to establish the meaning of these observations, we have 
grouped together the previously collected subsoil data from the zone under 
study. In this way a study of electrical soundings (GRGM-SGR Midi-Pyrenees- 
Mining Code n° 356-7-252) has in particular provided the resistivity and form 
of the basal surface of the surface limons. In addition, by grouping the obser- 
vations made on the wells of this zone by the BRGM-SGR Midi-Pyrenees we have 
reconstructed the pattern of the piezometric surface of the phreatic flow 
using observations made during 3 different periods (July, 1964; April, 1966; 
and May, 1966) (5,6). 

In figure 5 we have shown certain information relating to the obser- 
vations taken from IRC photography. 

Thus, if we compare the distribution of resistivities of the surface 
terrain with information obtained by balanced lighting, we generally note a 
good agreement betxj'een light zones and high resistivities. In figure 4, for 
example, we have sketched the development of the resistivity along two align- 
ments and the position of the dark zones meeting them. This thus leads us to 
interpret these local dark areas as being rich in clays and highly humidified. 
Their curvilinear 'concentric arrangement therefore outlines well the layout 
of former sandy meander cordons . 



Furthermore, the alignments of the areas of high surface resistivity 
and, simultaneously, relatively low limon thickness correspond well to those 
revealed under green light. We can interpret these wavy lines as former major 
beds of bodies of water at a stage of development prior to that leading to 
the sedimentation of the above mentioned cordons. In fact east of Beldour these 
cordons of average NNW-SSE alignment line up in some way with these long line 
patterns. 


Finally the extensive dark area obtained using red lighting is similar 
to the isopiezometric map. In fact, it seems to correspond to the area of low 
hydraulic gradient and thus of high permeability on the map. This is event more 
logical since it is noticeably lined up along the high points of the basal surface ; 
of the limons. Thus, in this zone, the vegetation, when there is any, would 
present a low rate of evapotranspiration compared with the very aerated state 
of the subsurface, while the bare topsoil shows the more sandy influence of the 
subsoil on its granulometry. 

These results thus make it possible for us to locate the theoretical 
outline that we reviewed for the formation of alluvial sediments : 

- An early stage with filling by more or less coarse materials 
corresponding, for example, to line A or B (figure 4) brought out by means 

of green lighting (it seems in effect that the filling of B may be more clayey 
because situated within the extensive dark zone) . 

- A second stage of sedimentation with the formation of very many 
point bars. The development of the channel as spared the first cordons formed 
and has rather "leaned" on them as is often evidenced by the parallelism between 
lines. The association of all these crescents creating a longline pile strongly, 
highly permeable characterized as a whole in red light. 

Conclusion : 

The localization of the position of former main beds within the 
alluvial deposits can be carried out using "false color" infra-red photographic 
documents. This treatment is improved by the use of a colored lighting system 
(red-green) with variable brightners coupled to a black and white TV circuit. 
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This preliminary study shows that it is possible to find not only 
the main axes of permeable formations but also the particular form of each • 
sandy point bar . 

This flexible rapid method therefore provides and inexpensive means 
of preparing a mechanical or geophysical prospecting study in an alluvial medium. 

N.B. We would like to thank Messieurs Angelier * Gourinard , Rey and Vincent 
of the Centre d' Etude et Amenagement des Ressources Naturelles (University 
Paul Sabatier) and Messieurs Aries and Flouzat of the Service de la Carte de 
la Vegetation (Centre National de Recherche Scientif ique) for the advice and 
material aid which they were willing to offer us. 
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Figure Captions : 


Figure 1 Diagram of the distribution of fluviatile deposits (The letters A, 
B, C, D, and E are referred to in the text). 

Figure 2 Effect of spectral shift on the appearance of the IRC photographs 
with or without colored lighting. 


Figure 3 Observation system for balanced lighting. C-Television camera - 
variable transformer. 

Figure 4 Main observations obtained by means of the system in figure 3 on 

photograph 74-IRC-2-5. Heavy broken line : alignments revealed under 
green light (A, B : see text) . Shaded area ; Extensive dark area 
under red light. Dotted area ; local dark areas under balanced lighting. 
The double line at top left and left marks the edge of the IRC photograph 

Figure 5 Main geophysical and hydrogeological information. Relative pattern of 
the piezometric surface of the phreatic layer. X and Y ; electrical 
sounding alignments showing in figures the amount of primary terrain 
resistivity. The black dots correspond to the dark areas obtained in 
figure 4. Heavy broken line : average alignments of shallow ground 
having higher resistivity. 



ABANDONNED 

BODY 

(DEAD ARM) 


FIGURE 1 



figure 2 



CONTRAST 


I’lGURT. 3 















FIGURE 4 


OF ^ Pm '^ pace is 

Qp:p;i5^' 


FIGURE 5 







CHAPTER 3 


A PRELIMINARY THEORITICAL STUDY IN ORDER TO ESTIMATE 
by REMOTE SENSING THE BIOMASS IN A POPLAR SITE 


G. FLOUZAT 


Abstract ; In this text, the theoritical approach to the study 
of the yield of poplar wood consists in determining the accuracy of the 
microdensitometric method of processing the aerial images of poplar sites. 

A reminder and a presentation of the active parameters show 
that it is the study of structures that should be undertaken rather than 
the study of the reflectance of the plantations, 

A thecaitical simulation model has been given along with an example 
of its application. Later, we will add to the report of the AGRESTE experiment 
the results of the thorought application of this model. 


This experimentation allows us to envisage a useful application 

of the data obtained. In an initial stage, the correlations with the ground 

biometric measures will indicate the efficiency of Remote sensing from a plane 
at a low altitude for the measurement of the workable part of biomass. 

At a second stage, the application of these results in a larger number of plots 

spreadover the whole site will increase the density of information of the 

ground -truth. 

The knowledge thus acquired of the region under study will serve us 
as a reference stick in judging the efficiency of the different systems when 
the resolution has diminished. A 1/20 000 map will thus be established and 
the data transmitted by the multispectral scanners (the aircraft Daedalus ; 
Landsat 2, MSS) will be compared. 


CESR : Centre d'Etude Spatiale des Rayonnement s 


TOULOUSE 


I. ACTIVE PARAMETERS IN THE STUDY OF POPLAR SITE 


The collation of parameters intervening in the measurement of 
the ratio of plantation cover gives the following list ; 


- Phenological stage of poplars 

- Spectral identification of crowns 

- Sun elevation on aerial survey time 

- Age of the plantation or yield class 

- Soil upkeep 

- plantation density 

- Plantation mode 

- Lopping of trees 

- Resolution of remote sensing imagery 

- Resolution of the analysis element of the imagery 

•- Angle between direction of tree lines - numerical scanning 
direction 

- Varieties cultivated 

- Soil water availability 

- Presence and fluctuations of water table 

- Chemical properties of soil. 


This accurate statement shows the factors which have to be known 
before one can look for acorrelation with the measures of biometric characte- 
ristics in the plots to be sampled. 

The influence of some of these parameters is eliminated or mini- 
mized by the conditions of the remote sensing process. This is true for the 
phenological stage, the position of the sun and the resolution of the image. 


In fact, the flight which GDTA has been comraissionned to do took place after 
total foliation has occured and at the time when the sun is at its highest in 
the region and at the apropriate date. The desired resolution corresponds in 
the films with the capacity of the photogramme trie camera in a photographic 
survey to the scale of l/lO 000. 

The parameter ’’the age of the plantation or the yield class” 
will be studied using the correlations obtained between the factors resulting 
from the biometric characteristics and the measurement of the cover rate 5 
which constitutes the basis of the model. 

We can then regroup the other active parameters into three 

factors : 

a) the conditions of plantation and upkeep 

b) the soil fertility 

c) the type of image processing 

In order to ascertain the influence of the first and third factor, 
a theoritical study was carried out before the aerial photographic data were 
obtained. 


II. SIMULATION OF THE MEASUREMENT OF COVER RATE 

It was particularly important therefore to knovj the accuracy of 
the measurement of the cover rate in the conditions selected for the image 
processing. This information constitutes in fact the only physical link with 
the biometric characteristics : 

- Diameter breast height 

- Basal area 

- Basis girth 

- Total height 

- Gutting height 

- Volume of wood cubic content 

For this, the theoritical model was built. It reproduces at 
a scale of 1/10 000 various plots where in most cases the crowns coincide 


with circles 



The figures on the next page show an enlargement (G = 2) of the model and the 

localization of the olots the characteristics of which are given in table I. 

allows ^ , . , , 

A reading of this, to affirm that the conditions of plantation and upkeep are 

simulated. 

This simulation was digitized with a step of 25 u by the micro- 
densitometer of the C.E.S.R, The availability of this numerical data obtained 
in 'teal dimension" allows one to study the role played by the image processing. 
Xn this report an example of an attempt to establish the accuracy of the measu- 
rement of cover rate is given for plot n° 1. A prerequisite of this analysis 
was the definition and the application of the process represented by the diagram 
(Software of C.E.S.R. data processing group). Two examples of intermediary 
stages are given in the figures "Display of plots n® 1, 5 and 6" and 

\ 

"Histogram of plot n° 1. 




Theoritical model used 


(Simulated scale ; 1/5 000) 











TABLE I 


PLOT CILARACTERISTICS OF THE THEORITICAL MODEL 


Plot number Diameter 


Crown 

girth 


15,70 m 


15,70 m 


Crown Plantation 


19,63 m2 7 X 7 m 


Density/ha Cover rate Angle to Observations 
ratio 


15,70 m 19,63 m2 7 x 7 m 


15,70 m 19,63 m2 7 x 7 m 


19,63 m2 7 X 7 m 


15,70 m 19,63 m2 6 x 6 m 


13,43 m 14,38 m2 6 x 6 m 


3,14 m 0,78 m2 j 7 x 7 m 


6,28 m 3, 14 m2 7 x 7 m 


10,99 m 9,61 m2 7 x 7 m 


21,19 m 35,76 m2 7 x 7 m 


21,19 m 35,76 m2 7 x 7 m 


15,70 m 19,63 m2 6 x 8 m 


15,70 m 19,63 m2 6 x 8 m 


15,70 m 19,63 m2 6 x 8 m 


15,70 m 19,63. m2 6 x 8 m 


25, 12 m 50,24 m2 7 x 7 m 






















































































A series of element analysis taken on plot n° 1 and carried out 
in accordance with the procedure of the diagram is expressed on .the table II 


TABLE II - ANALYSIS ON PLOT N° 1 
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Display! of theoritical plots, n° 1, 5 and 6 



Histogram of plot n° 1 





DIAGEIAM OF MODEL ANALYSIS 



The distribution of these measuretnents of cover rate are repre- 
sented in the fig. "Distribution of the measurements of the cover rate". 

This simulation presumes that the parameter "spectral discrimination of the 
crovms" is really known. This calls for a special study on significant 
pass-bands in this kind of remote sensing. 










Distribution of the measurements of cover rate 


Number of 
measurements 



^ 7c cover 


Ml A0,1 
M2 = 41,1 


Cover rate of the model 
Average of cover rate measured 


The difference between Ml and M2 is satisfactory and allows one 
to presume that measurements on aircraft data will be accurate enough to be 
applied in the model o 

However, one notices a general shift in the measurements in com- 
parison with the real value. This is due to the size of analysis element. 

The sides of the squares or rectangles in which the histogram is effected are 
slightly inferior to the distance betv;een two lines of plantations. Thus, 
this increases the ratio givingthe cover rate. In order to obtain values 
which are closer the reality, one must use surface elements the sides of which 
are exactly proportional to the spacing of the plants. 



This theoritical study will thus allow an appreciation of the 
influence of the parameters which are directly linked to the processing of 
data obtained by remote sensing as regard the study of poplar sites. It has 
the added interest of revealing all the factors which play a part in this 
attempt to simulation. 

This has been made easier by undertaking a first inventory of 
poplar plantations all over the middle Garonne area (about 50 Icms from 
Toulouse to Castelsarrazin) , This map is at the scale of l/20ooo and obtained 
by photointerpretation from I.G.N, aerial photographical missions of 1969 
and 1971 (panchromatic emulsions). This document was elaborated by including 
subdivisionsin the plots where the growth of the plants is not homogeneous. 


This work will enable us to tackle the problems of biomass 


estimation in favorable conditions. 


CHAPTER 4 


PREMIERE PARTIE DE lA CARTOGRAPHIE 
DES PEUPLERAIES DE LA MOYENNE VALLEE DE LA GARONNE 
(Test-site n^S du Programme AGRESTE) 


G. FLOUZAT 


Resume : 

La premiere partie de la carte des peupleraies de la moyenne 
vallee de la Garonne est px'esentee a son echelle originale 
(E = 1/20 000). E' ouite, le programme d 'utilisation de cette 
carte et des photographies infra-rouge couleur est decrit, 

Les phases successives doivent amener une connaissance du site 
directement utilisable et une reference pour le traitement des 
donnees acquises par les autres moyens de teledetection sur le 
mime site. 


I - CARTOGRAPHIE DES PEUPLIERS 


La carte des peupliers a ete realisee par deux photointerpreta- 
tions successives ; 

- en 1974, la premiere utilisait des photographies panchromatiques 
(echelle : 1/25 000) des missions de I'Institut Geographlque 
National (I, G.N.) (Montauban 1969 j Toulouse-Ouest 1971), 

- en Decembre 1975, la seconde utilisait des leur disponibilite 
les photographies infra-rouge couleur (echelle l/lO 000) du 
Groupement pour le developpement de la Teledetection Aerospa- 
tiale (G.D.T.A.) (Mission "peupliers de la Garonne" du 

01 - 05 - 1975). 


GESR : Centre d'Etude Spatlale des Rayonnements 


TOULOUSE 



Cette interpretation touche uniquement les peupliers et ne repre- 
sente pas les autres formations boisees. La zone cartographiee couvre 
la vallee de la Garonne sur une vingtaine de kilometres, Les quatre sec- 
tions de carte necessaires correspondent au tiers de 1’ ensemble du test- 
site n° 5 du programme AGRESTE, 

Le critere de photo interpretation est celui defini par 1‘lnstituto 
per la Sperimentazione alia Pioppicultura" avec quatre classes de recou- 
vrement : 


Glasse 

1 

; 0 

a 

5 ' 

Glasse 

2 

; 5 

a 

25 

Glasse 

3 

; 25 

K 

a 

75 

Glasse 

4 

; + 

de 

75 


Get inventaire repertorie 451 parcelles de peupliers dont la 
repartition est donnee en II - 1. 

Une observation a ete faite en comparant les photointerpretations 
des missions 1969 - 1971 et 1975 ; elle fait apparaltre que norabre de 
peupleraies exploitees depuis 1971 n’ont pas ete replantees en 1975. 

Les pages suivantes presentent la position des quatre sections 
de carte par rapport a 1 ’ensemble du test-site j puis chacune de ces sec- 
tions a I’echelle originale (1/20 000). 




LOCALISATION DU SITE N« 5 



DU PROGRAMME AGRESTE 

ET DE LA PREMIERE PARTIE DE LA CARTOGRAPHIE 
DES PEUPLERAIES 


Le trait double delimite la region etudiee (Moyenne Vallee 
de la Garonne) et les rectangles correspondents aux quatre 
sections cartographiees. 




SECTION 1 


LEGSNDE 


Pcupleraies 


Recouvrecent 
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SECTION 2 


LEGENDS 


Pcupleraies 


Rccouvrecent 
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5 % 
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a 

25 % 
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de 

75 % 









SECTION 3 




LEGEKDE 
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SECTION 4 


Peupleraics 
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II - PROGRAMME D 'UTILISATION DE LA CARTE ET DES DONNES IRC 


L 'acquisition de cette verite-terrain a deux buts ; 

a) - connaltre I'etat actuel de la oiture du peuplier dans cette 

region 

b) - disposer d'une reference permettant une comparaison avec les 

informations apportees par d'autres capteurs de teledetection 
et ainslj connaltre leurs limites respectives dans ce type 
d' experiment at ion. 


1. Resultat immedlat et premiere conclusion 

La repartition des 541 parcelles repertoriees est donnee 
dans le tableau suivant ; 


Nombre de parcelles 


classes 

Sections 

1 

2 

3 


1 

0 

1 

3 

14 

2 

10 

14 

17 

27 

3 

22 

24 

50 

97 

4 

13 

15 

37 

87 


Ce ccmptage sera suivi dans une prochaine phase de la mesure 
des surfaces occupees par cheque classe de recouvremetit . 

Pour toutes les parcelles, le resultat est le suivant ; 


classes 

1 

2 

3 


Nombre de parcelles 

55 

64 

107 

225 

Taux 

127o 

147. 

237. 

507. 


Ce tableau entraine une premiere conclusion relative a la 
P^^'^^fsion globale de production de bois de peuplier dans cette region. 



























Bien que la mesure des surfaces coorespondantes ne solt pas encore faite 
definitiveraent, il semble qu'll y ait une baisse tres nette de la planta 
tion du peuplier au cours des dernieres annees. Il est done probable 
qu'il y aura une baiss e de la production de bois de peuplier dans cette 
zone lorsque la classe 2 et la classe 1 actuelles seront en age d'etre 
exploitees. 


2. Echantillonnage et mesures de terrain. 

Ces operations ont pour but d'apporter les informations 
necessaires a 1 ' etablissement de relations entre la structure du "toit" 
des peupleraies et leur volume de bois. Elies s'appuient done sur une 
conc eption qui ne fait pas appel a la mesure d*une reflectance. 
L'utilisation des photographies infra-rouge couleur est justifiee dans 
ce cas par la resolution spatiale qu'elles fournissent. 

Ainsi, les caracteristiques de cette detection, simulees 
theor iquement montrent que la technique d' etude du taux de couverture 
donne une mesure a + 3 % pres. Or, la periode significative de I'accrois- 
sement du recouvrement se situe entre 5 et 80 °L, 

Ces deux conditions de travail font que la methode d'analys^ 
microdensitometr ique des peupliers permet de distinguer un nombre de 
classes suffisant pour la correlation avec les stades de croissance 
interessants . 

Ces stades seront echantillonnes selon les criteres suivants 

- homogepeite des parcelles 

- representativite et complementaf ite en age des plantations 

- distance au point central de la photographie servant a 
1' analyse n' introduisant pas de distorsion. 


Les leves biometriques seront effectuees sur les peupleraies 
retenues. La technique utilisee sera la photographie hemispherique ou 
huit mesures sont signif icatives dans chaque prise de vue. 
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Arbre me sure 
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Appareil photo graph ique 


Disposition des 8 arbres etudies 


Pour un arbre, les mesures suivantes sont effcctuees sur la 

photographie : 

- surface terriere (avec une part de contr61e au sol) 

- hauteur a la decoupe 

- cubage par billons 

A cela s'ajoute la mesure au sol de la surf ace basale . 

Les plantations etant homogenes, trois prises de vues suffiront 
pour definirle volume moyen des arbres de la plantation etudiee. (Soient 
24 mesures par parcelle). 

Ce procede implique un plan de charge operationnel fixant 
le nombre de parcelles a retenir vers 50. Cinquante mesures moyennes seront 
done tirees de 1 200 cubages. Get echantillonnage affectera ainsi 11 7. 
des peupleraies de cette partie de la carte. 


3* Exploitation 

Ces donnees levees sur le terrain seront ensuite mises en 
relation avec les mesures des taux de couverture des memes parcelles par 
microdensitometrie . Cette mesure fournira la surface du houppier vue a la 
verticale. Cette surfam est dependante des conditions de plantation. 

Trois cas principaux sont representes dans la moyenne vallee de la Garonne : 


- Plantation "carree" : 7m x 7m (204 plants/ha) 

- Plantation "carree" : 6m x 6m (278 plants/ha) 

- Plantation "rectangulaire" : 6m x 8m (208 plants/ha) 


Les trois graphiques suivants donnent les relations taux de 
couverture “Surf ace du houppier vu a la verticale* 

Relation Taux de couverture - Surface du houppier pour une plantation 7 x 7m 

Surface du houppier 
Sh m®A 



Relation Taux de couverture - Surface du houppier pour une plantation 6 X 6m 
Surface du houppier 



5 25 50 75 



Relation Taux de couverture - Surface du houppier pour une plantation 6 x 8ra : 


Surface du houppier 



L'obtention de cette mesure du houppier servira a la recherche 
des lois reliant ; 

Taux de couverture a Surface terriere 
Taux de couverture a Voluir ^ hois 

Les mesures au sol s'ajoutant aux correlations eventuelles 
entre mesures de teledetection et mesures de terrain, deux types de 
resultats sont a esperer ; 

Une table Surface terriere / Volume de bois 

Une table Taux de couverture / Surface terriere 

11 va de soi que 1 ' echantillonnage de base est conqu dans 
des conditions ecologiques homogenes pour eliminer les perturbations dues 
a des croissances de physionomies differentes. 

A 1 ' etablissement des deux tables, s'ajoutera la relation 
entre le volume de bois exploitable et I'age de la plantation puisque 
cette information sera collectee au cours de 1 'etude des parcelles echan- 
tillonnees. Ctci devait conduire aux deux coiinaissances suivantes ; 


- evolution de 1 ’accroissement annuel d*une peupleraie 

(I’Sge donnant la function du temps) 

- prevision de. la production disponible d»une z6ne 

(la carte fournissant la mesure de surface) 


Dans le cadre de ces operations d’utilisation de la Carte 
des Peupleraies et des donnees IRC, des echantillonnages diffdrents 
entrepris ultdrieurement seraient souhaitables pour juger de la variation 
des tables en function des conditions pedologiques d’ installation des 

plantations. 



CHAPTER 5 



COMPARAISON ET PREMIERE INTERPRETATION 
DE DIFFERENTES DONNES DE TELEDETECTION 
SUR UNE MEME ZONE DE PEUPLERAIES 


G. FLOUZAT (CESR) et J. DAGNAC (SCV) 


Resume : Un bilan des donnees disponibles pour la realisation de 

I'objectif "Peupliers" du programme Agreste dans la Vallee 
de la Garoune montre que 1 ’ information utilisable efficace- 
ment est contenu dans un petit nombre de scenes. 

La date tardive de reception de ces donnees n’interdit pas 
de faire une comparaison entre les images obtenues aux 
differents niveaux de teledetection. L ' interpretation four- 
nit alors une recapitulation des renseignement s que l‘on 
peut obtenir de chaque type de document. 


I - DONNES DISPONIBLES ET DONNES COMPAREES 


Les donnees disponibles pour 1 ' etude du site n° 5 
du programme Agreste sont indiquees dans le tableau de la page survante, 
Elies representent la quasi-totalite de 1 • information et les differentes 
conditions d’ acquisition. Pourtant, le volume de traitement necessaire 
dans cet objectif est limite. En effet, deux scenes Landsat 2 d'excel- 
lento qualite representent toute I’imagerie interessante a petite cchelle 

CESR : Centre d’Etude Spatiale des Rayonnemenn s 

TOULOUSE 

SCV : Service de la Carte de la Vegetation 
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TABLEAJ DES DONNEES DISPONIBLES 


V ehicules 
capteurs 


Daedalus MSS 


IRC Photo. 


Altitude 


Bandes 


Landsat 2 MSS 02-06-75 910 Kin 



01-05 


1 500 m 0j55-0,80 


Resolution Scenes 

du Site n®5 


1/3 Nord 


2/3 Sud 


1/4 Sud 


" 2/3 Sud 


2/3 Sud 


1/2 Nord 


3/4 Sud 


Qualite Observations 


Mauvaise 
Nuages 



Moyenne 
Qq nuages 


Nuage-Site Bande magnetique 
utilisable demandee 



Bonne 


Bonne 


Bande magnetique 
demandee 


Site entier Bonne 


1/10 Sud 


Nuage-Site 

utilisable 



Site entier Bonne 


Site entier Bonne 


Bandes numeriques 
possibles sur les 

zones interessantes 


_ j 1 n 49 • 0 42-0.45 i 0,45-0,50 J 0,50- 0)55 j 0)55-0)60 

* Bandes passantes du Scanner Daedalus : 0)3c-0) > > > 

0)60-0.65 ; 0)65-0)70 ; 0)70-0,80 ; 0,80-0,90 ; 0,90-1,10 





















































Les autres donnets (imagerie a grande echelle ; MSS Daedalus, 
Photographie IRC) dont le trattement est necessaire sont egaleraent en 
quantite ralsonnable puisqu'un echant illonnage regional eliminera les 
zones redondantes. 

Une partie cartographiee du site n° 5 a ete choisie 
pour presenter un exsnple de comparaison d’imgges a des altitudes tres 
dissemblables (sections 3 et 4 de la carte au 1/20 000), 

Les documents retenus sont les suivants : 

- Carte des peupliers dans la partie Sud n® 5 

(reduite a 1' echelle : 1/100 000) 

- Reduction d'un assemblage noir et blanc des photographies 

infra-rouge couleur. (Echelle : l/lOO 000) 

- Reduction d'une visualisation de donnees Daedalus. 

(Echelle : l/lOO 000) 

- Agrandissement de la region cartographiee avec une scene 

Landsat 2. (Echelle ; l/lOO 000) 

L’utilisation dans une premiere interpretation et la 
preparation de ces montages aboutissent a la description des apports 
respectifs des principaux types de teletctection. (Les bandes passantes 
choisies sont equivalences). 


II - APPORTS RESPECTIFS DES TROIS CAPTEURS 


1. Photographies infra-rouge couleur 

La mission aerophotographique au l/lO 000 realisce 
avait pour but de fournir un inventaire detaille des peupleraies. Les 
caracteristiques de resiution des images obtenues ont done donne lieu 
a une cartographic des surfaces occupees par quatre classes d'age. Les 
rapports entre les differences classes d’gge pour etablir une prevision 
et les conditions d' execution en sont fournies au chapitre precedent : 
"Premiere partie de la cartographie des Peupleraies de la moyenno vallee 
de la Garonne". Une reduction de cette carte au 1/100 000 est presentee 
comme reference pour aider a la comparaison globale des images des diffc- 


rents capteurs. 











Les photographies infra-rouge couleur apportent d'autres 

renseignement s : 

- Caracteristiques du type de plantations 

- Extension de la zone de faible distorsion 

- Etat et entretien du sol 

- Etat de feuillaison 

- Taux de couverture. 

REDUCTION D'UN ASSEMBLAGE NOIR ET BLANC DES PHOTOGRAPHIES 
INFRA-ROUGE COULEUR SUR LES SECTIONS 3 ET 4 



- Caracteristiques du type de plantation. 

Les me£,ures photogrammetriques permettent, apres le calcul 
precis de 1‘echellej de caracteriser la peupleraie par : 

le mode de plantation : carre ou rectangulaire 
les distance entre les lignes de plantation 
' la denslte de pla.itation. 

Ainsi, le calcul du taux de couverture peut §tre fait en 
tenant compte de ces parametres pour exprimer la surface du houppier. 


- Extension de la zone de faible distorsion. 

I’examen stereoscopique des photographies indique jusqu'a 
quelle distance du nadir la distorsion ne risque pas d’entrainer une erreur 
dans la mesure du taux de couverture. En effet, lorsque cette distance aug- 
mente, la mesure de la surface des houppiers est faussee par 1 * introduction 
dans le recouvrement de point appartenant a des peupliers vus "de profil" 
ou en perspective. La zone de faible distorsion doit done indiquer les lieux 
ou les arbres sont photographies "verticalement". 


- Etat et entretien du sol. 

La presence ou 1’ absence d'un tapis herbace© dans une 
plantation de peupliers revet une grande importance dans les peupleraies 
jeunes. C’est dans ces classes que la confusion ei.tre la couronne de I'arbrc 
et la vegetation au sol risque de nuire le plus a la determination du taux 
de couverture. Dans des plantations plus ageea, le recouvrement devient plus 
grand et 1' erreur sur sa mesure dirainue en meme temps que la faible partie 
de sol visible a la verticale se reduit et se trouve masqiieepar les ombres 
des peupliers. 

Les photographies infra-rouge couleur indiquent nettement 
I'etat du sol. Il est done possible de connaltre les parcelles ou les 
les mesures des houppiers devront ctre controlees. 



- Etat de feuillaison. 

La date de la mission aerophotographique est le ler mai 
'1975. Dans les conditions biocliinatiques du debut de cette annee 1975, 
les peupliers n*ont pas toujours une feuillaison complete. Il est ainsi 
possible de distinguer 3 stades d*avancement de cette phase phenologique : 
debourrage 
feuillaison moyenne 
feuillaison complete. 

Il apparait generalement sur les photographies que ces 
differents stades ne sont pas presents dans une merae parcelle, Dans quel- 
ques cas, il a ete constate qu'ils permettaient de distinguer deux variete 
tres proches dans I'espace. Une etude est done entreprise pour verifier 
ou detruire I’hypothese que ces stades d'avancement de la feuillaison cor- 
respondent a des differences varietales ou clonales des peupliers. 


- Taux de couverture. 

L'echelle des photographies (1/10 000) a ete definie en 
function des enseignements tires de 1 ‘analyse de la maquette theorique dont 
un example est donne dans ce rapport au chapitre "A preliminary th eoriti- 
cal study in order to estimate by remote sensing the biomass in a poplar 


site". 


L’ analyse microdensitoiiietrique du recouvrement sera done 
appliquee aux peupleraies dans des conditions methodologiques connucs. 

De plus, la definition des caracter istiques du type de plantation foiirnit 
les parametres necessaires pour indiquer la surface du houppier en fonction 
du taux de couverture. Enfin, il faut noter que la risolution spatiale du 
film infra-rouge couleur aerien utilise est voisin du pas de digitalisation 

de 1 ' image ( 25 cm) . 


2. Radiometre a balayage Daedalus 


Une couverture aerienne a basse altitude sur la moyenne 
vallee de la Garonne a ete realisee le ler mai 1975 par le radioraetre a 
balayage Daedalus puisque celui-ci etait embarque dans le meme a^fion que 


omswM, page s 
POOR 



la camc'ra Willi cliarg.'e avoc Ip film infra-rouge coulcur. La comparaison 
des deux reductions au 1/100 000 des assemblages IRC et Daedalus montre 
une resolution spatiale plus taihlc du radioraetre (environ 3 m) . 


reduction D’UNE visualisation DE DONNEES DAEDALUS 
SUR LES sections 3 ET A 


Echelle : 1/100 000 



Canal 7 : 0,f5 - 0,70 u 


Canal 10 : 0,90 - 1,10 u 




Par ailleurSj iL serable que plusieurs canaux combines ne 
permettent pas forcement de discriminer les peupleraies, L‘ ident if icaCion 
par la reconnaissance de structure est aleatoire du fait de la valeur du 
rapport objet/resolution spatiale. La discrimination raultispectrale ne leve 
pas toujours I'ambiguite par rapport aux formations boisees appartenant au 
cortege floristique des associations hygrophiles et les autres plantations 
(fruitieres) sont egalement difficiles a separer. 

Dans les jeunes peupleraies, 1* importance de I’etat et 
de I’entretien du sol est primordiale. Lorsque le sol est nu, la resolution : 
est telle que la petite taille du houppier ne modifie pratiquement pas la ' 
reponse qu * enregistre le radiometre. Lorsque le sol est convert d'un tapis 
herbace, le radiometre reqoit une energie correspondent a un toit continu 
de vegetation. Dans les deux cas, les classes d'age jeunes sont affectees 
par erreur, soit au sol nu, soit a des classes plus agees. 

Il faut done s'attendre a des difficultes dans I'utilisatic 
des classifications multispectrales pour la discrimination des peupleraies. 

Dependant, les donnees numcriques du radionietre a balayage 
Daedalus seront du plus grand interet lorsqu'il deviendra interessant 
d’etudier la perte d* information en fonction de la degradation de la reso- 
lution spatiale. Cette operation, qu'il est prevu de mener conj ointement 
avec une analyse identique sur la maquette theorique, pourrait apporter 
d’utiles renseignements sur la signification de 1 ' inf ormation contenue dans 
les donnees LANDSAT 2 dont la resolution est de 80 m« 

3. Radiometre a balayage LANSAT 2 


Les donnees du scanner multispectral du satellite LANDSAT Ti 
fournissent des informations a petite echelle, conforraes a ce que laissnit 
suppose!" la mission Daedalus. 

Ainsi, la discrimination ne semble pas pouvoir separer les ; 
diff e-rentes formations boisees (Peupliers, Fruitiers, bois du bord des eaux)i 
Par centre, 1 ’ identification des cultures et de 1 'hydrographle est efficacc.: 
Un renseignement interessant est obtenu dans cette iinageric 
grSce a la vue synthetique que procure 1 ' observat ion a haute altitude et la 
faible resolution : il s’agit de la delimitation par photointerpretation 


AGRANUISSEMENTS DE LA REGION CARTOGRAPHIEE 
DU SITE " PEUPLIERS - GARONNE " 


(Scene LANDSAT 2 du 07 - 07 - 1975 ; E - 1/200 000) 








de la 2 one proche de la Garonne contenant la quasi-total ite des formations 
boisees de la vallle. (Un example de cette cartographie est donne sur la 
partie cartographiee du site.) L'interet de cette infomation est qu'elle 
correspond vraisemb lab lament a la basse plaine. Cette zone est done celle 
qui reunit les conditions ecologiiues optimalespur le peuplier. 


CARTOGRAPHIE DE LA ZONE PROCHE DE LA GARONNE 
CONTENANT LES FORMATIONS BOISEES DANS LES SECTIONS 3 ET 4 



Cette donnee apporte ainsi la possibilite de mesurer le 
potentiel de production de la region par les peupleraies, En effet, cette 
zone est fomiee de differents elements ; 

Peupleraie s 

Formations boisees naturellcs 
Fruit iers 

Cultures et Prairies 

Surfaces en eau (Garonne et gravieres) 




La surface totale de la zone est facile a connaltre ainsi 
que 1 'ensemble des formations boisees et des surfaces en eau. La difference 
permet de co maitre 1 ' importance des cultures et des prairies. 

Par la cartographic, les terres occupees par les peuple- 
raies sont connues et indlquent la proportion qu'elles representent par 
rapport a 1 'ensemble de la roue). 

Deux types de potent iels de production seront ainsi 

definis ; 

1. La quantite de bois de peuplier que fournitait ce 
milieu si les plantations occupaient 1' ensemble des 
surfaces boisees. 

2. La quantite de bois de peuplier que fournlrait ce 
milieu si les plantations occupaient les surfaces 
actuellement en peuplier et les surfaces en boisemcril 
naturel . 

(L'hypoUliese d'oecupation des terres prcalablc a la 
definition de ce potentiel exclut 1 ' installation du peuplier sur les 
terres utilisees actuellement par 1 ' agriculture . ) 


En conclusion, Ic prmier exarnen des differences donnees sur unc 
mSme zone de peupleraies fait apparaitre les posslblilites et les liiuites de 
chaque type de donnees. Ces resultats utilisables sont lies aux conditions 
de leur acquisition et montrent une certaine complemcntar ite . 



